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I.  INTRODUCTION  AND  BASIC  EQUATIONS 

Adiabatic  shear  banding  is  a  localization  phenomenon  that  occurs  dur¬ 
ing  high  rate  plastic  deformation  in  many  materials.  Since  the  heat  gener¬ 
ated  by  plastic  flow  usually  tends  to  lower  the  flow  stress,  the  response 
curve  for  adiabatic  homogeneous  deformation  will  fall  below  the  isothermal 
response.  If  this  thermal  softening  is  stronger  than  both  work  and  rate 
hardening  together,  then  eventually  a  peak  stress  occurs  at  a  critical 
value  of  strain,  followed  by  decreasing  stress  for  further  increments  of 
strain.  Such  a  situation  is  generally  unstable  and  presents  an  opportunity 
for  strain  localization  to  occur.  In  recent  years  there  has  been  consider¬ 
able  interest  in  developing  a  quantitative,  dynamical  theory  of  adiabatic 
shear  banding.  Much  of  the  recent  literature,  as  well  as  some  of  the 
pioneering  works  have  been  listed  and  discussed  briefly  by  Clifton,  et  al- 
In  this  paper  we  present  further  calculations  of  the  kind  given  previously 
by  Wright  and  Batra.2  in  particular,  we  show  the  effect  of  changing  the 
constitutive  rate  response  and  the  magnitude  of  the  perturbation.  We  also 
examine  the  effect  of  adding  a  gradient  or  dipolar  response  to  the  consti¬ 
tutive  equations,  at  least  for  early  times. 

In  order  to  concentrate  on  fundamentals,  consider  one-dimensional 
shearing  of  a  block  of  material  that  lies  between  the  boundaries  y  =  ±h. 

The  deformation  is  assumed  to  be  given  completely  by  horizontal  shearing. 
Thus  the  velocity  field  may  be  written 

x  =  v(y,t) ;  y  =  Q;  z  =  0.  0) 

As  used  in  Reference  2,  the  balance  and  constitutive  laws  for  this 
case  are  given  here  in  nondimensional  form. 


Momentum: 


Energy : 


Elastic  Response: 


s’y  =  pv 


e  =  *0,yy  +  syp 


s  =  u(v,y  -  Yp) 


Reference  Plastic  Response:  tc  =  (1  +  )n 

^0 


Work  Hardening; 


Yield  Surface: 


Wp  *  K*  =  SV 


[ s 1  =  (1  -  ae)(l  +  b 1 Yp|)  < 


Boundary  conditions  are  v(±1,t)  =  ±1,  6,  (+1,t)  =  0.  In  these  equations 

the  shear  stress  is  s,  ma3s  density  is  p,  particle  velocity  is  v,  tempera¬ 
ture  rise  is  0,  thermal  conductivity  is  k,  plastic  strain  rate  is  Y  ,  shear 

modulus  is  p,  work  hardening  parameter  is  x,  and  plastic  strain  in  a  refer¬ 
ence  test  is  ,  The  comma  denotes  partial  differentiation  with  respect  to 


rt 
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the  spatial  coordinate  y,  and  the  superimposed  dot  denotes  partial  differ¬ 
entiation  with  respect  to  time  t.  In  the  energy  equation  plastic  working 
acts  as  a  source  term.  A  linear  elastic  response  is  assumed  as  well  as  the 
additive  decomposition  of  shear  strain  into  elastic  and  plastic  parts.  The 
reference  plastic  response  is  intended  to  be  simply  a  curve  fit  to  a  slow 
isothermal  test.  The  work  hardening  parameter  is  assumed  to  depend  only  on 
the  plastic  work  no  matter  what  the  rate  of  the  test.  Finally,  the  yield 
surface  is  taken  to  depend  on  plastic  strain  rate,  as  well  as  stress  and 
temperature.  When  (2)g  is  solved  for  Yp  as  a  function  of  s,  0 ,  and  <,  this 

last  condition  is  seen  to  be  simply  an  overstress  rule  for  plastic  flow 
when  the  stress-temperature  point  lies  outside  the  static  yield  surface. 

The  form  chosen  here  is  similar  to  one  due  to  Litonski.  3  Note  that  the 
thermal  softening  depends  only  linearly  on  temperature. 


In  Reference  2  the  nondimensional  variables  and  parameters  are  related 
to  their  dimensional  (barred)  counterparts  as  follows: 


y  =  y/h 

t  =  tv0  V  =  v/(hY0)  S  =  1/kq 

o 

V 

> 

o 

1  Q. 

1  CD 

II 

CD 

A 

II 

A  1 

A 

O 

2  2 

y  *  y  =  $  p  a  (ph  y0  )/<0 

k  =  k/(pcvYQh2) 

p  »  »/<0 

a  =  (a< 0)/(pcv)  b  =  Byq 

V  =  V  ^0 

where  v 
'o 

=  v(h,t)/h  is  the  average  applied 

strain  rate  between  y  =  ±h,  < 

is  the  initial  yield  stress  in  the  reference  test,  cy  is  the  specific  heat 
at  constant  volume,  and  a,  b,  tpQ,  ra,  and  n  are  material  constants. 


II.  HOMOGENEOUS  SOLUTIONS  AND  PERTURBATIONS 

If  we  set  v  =  y  and  assume  that  s,  6,  and  <  are  independent  of  y,  then 
equations  (2)  reduce  to  a  set  of  ordinary  differential  equations  in  time 
with  s,  e,  and  <  as  dependent  variables.  For  initial  conditions  we  assume 

s(0)  =  1,  0(0)  =  0,  and  <(0)  =  1 ,  so  that  yielding  begins  at  the  initial 

time.  For  the  assumed  visco/plastic  flow  law,  the  stress/strain  response 
is  a  smooth  curve,  which  rises  at  the  elastic  slope  initially  and  then,  as 
plastic  flow  increases  and  the  temperature  rises,  the  curve  bends  over, 
passes  through  a  single  maximum,  and  finally  decreases  with  further 
increments  of  strain.  This  type  of  calculated  behavior  is  well  known  in 
the  literature,  e.g.  see  References  4,  5,  6,  and  7.  Figure  1  shows  typical 
results.  For  curve  A,  the  nondimensional  parameters  are 

p  =  3.928  x  10-5,  k  =  3.978  x  10-3,  a  =  0.4973,  v  -  240.3 

n  =  0.09,  <h  =  0.017,  b  =  5  x  106,  m  =  0.025 

5 

and  for  curve  B  all  the  numbers  are  the  same  except  b  :  5  x  10  and 
m  =  0.02.  Also  shown  in  the  figure  are  the  reference  stress/strain  curve 
(a  =  0,  b  =  0)  and  the  isothermal  response  (a  =  0)  corresponding  to  A  (case 
B  is  similar) . 
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Perturbations  to  the  homogeneous  response  have  been  introduced  in  the 
following  manner.  Just  prior  to  the  occurrence  of  peak  stress,  the  temper¬ 
ature  was  modified  by  adding  a  symmetric  bump  at  the  center.  Two  cases 
have  been  considered  as  shown  in  Figure  2,  the  larger  perturbation  being 
five  times  higher  and  broader  than  the  smaller  one,  but  similar  in  shape. 
With  the  new  temperature  distribution  given,  the  stress  was  recalculated  so 
that  (2)g  is  still  satisfied,  all  other  variables  being  held  fixed.  Then 

the  full  problem  using  all  of  (2)  and  the  stated  boundary  conditions  was 
restarted  with  the  new  field  variables  as  initial  data.  After  casting  the 
equations  into  weak  form,  calculations  were  made  using  the  finite  element 
method,  e.g.  see  Becker,  et  al  .3 

Case  B  with  the  smaller  temperature  bump  as  a  perturbation  was  previ¬ 
ously  described  in  some  detail  in  Reference  2.  Figures  3  through  7  illus¬ 
trate  and  compare  the  main  features  of  the  response  for  all  cases  consid¬ 
ered  to  date.  Because  of  the  constant  velocity  maintained  at  the  bound¬ 
aries,  average  strain  rate  in  the  scaled  variables  for  the  strip  y  is 

ave 

exactly  equal  to  one.  This  fact  allows  direct  comparisons  between  the 
homogeneous  and  inhomogeneous  cases.  Figure  3  shows  the  stress  at  the 
center  of  the  strip  as  a  function  of  average  strain  (or  time).  The  pertur¬ 
bation  was  added  at  the  point  marked  I,  and  the  peak  of  the  homogeneous 
response  is  marked  P.  For  the  small  temperature  bump  the  stress  follows 
the  homogeneous  response  until  well  past  the  point  P,  and  then  drops 
rapidly  as  the  shear  band  accelerates  and  localizes.  The  stress  at  other 
points  of  the  strip  is  nearly  the  same,  but  lags  slightly  in  time.  Note 

that  in  real  time  (for  yq  =  500  s  ,  say)  the  stress  collapse  is  delayed 

approximately  125  us  past  peak  stress  P,  and  the  rapid  acceleration  itself 
takes  more  than  20  ps.  This  behavior  is  entirely  similar  for  both  cases  A 
and  B.  The  case  of  the  larger  temperature  perturbation  contrasts  sharply. 
In  this  case  the  stress  collapse  begins  even  before  point  P. 

Figures  and  5  show  the  evolution  of  temperature  and  plastic  strain 
rate  at  the  center  and  edge  of  the  band  in  comparison  to  the  homogeneous 
case.  As  is  to  be  expected,  the  behavior  parallels  the  stress  response. 

At  first  the  temperature  and  strain  rate  deviate  3lowly  from  the  homogene¬ 
ous  response,  but  when  the  stress  collapses,  they  both  rise  steeply  in  the 
center.  At  the  edge,  the  temperature  levels  off  to  a  plateau  while  the 
plastic  strain  rate  (not  shown  for  the  edge  position)  drops  towards  zero. 

At  first  it  appears  that  the  temperature  and  plastic  strain  rate  change  so 
as  to  compensate  each  other  without  appreciably  affecting  the  stress 
response,  but  eventually  thermal  softening  wins  out  in  the  center.  At  the 
edge  the  stress  drops  due  to  momentum  transfer  from  the  center  and  effec¬ 
tively  quenches  both  the  temperature  rise  and  the  plastic  deformation. 

Figures  6  and  7  show  cross  sections  of  temperature  and  plastic  strain 
rate  for  case  A  (with  the  larger  temperature  perturbation)  at  several  times 
during  the  rapid  transition.  The  most  notable  feature  here  is  the  narrow¬ 
ing  of  the  most  active  region  of  deformation.  This  is  particularly  evident 
in  the  cross  sections  of  plastic  strain  rate.  Note  that  for  the  latest 
time  the  plastic  strain  rate  in  the  center  of  the  band  is  nearly  80  times 
the  average  applied  strain  rate. 
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since  y  =1).  Homogeneou 
ave 


Temperature  rise  for  case  A.  Upper  curves  are  for  the  center  of  the  band 
lower  curves  for  the  edge.  Solid  curve  is  the  homogeneous  case. 
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III.  MODIFICATIONS  FOR 


As  a  shear  band  forms,  it  is  evident  that  very  3teep  local  strain 
gradients  must  develop.  Therefore,  it  has  seemed  worthwhile  to  reformulate 
the  problem  in  terms  of  a  dipolar  theory  of  plasticity.  This  has  been 
accomplished  in  a  straight  forward  manner  by  modifying  the  theory  due  to 
Green,  Mclnnis,  and  Naghdi9  to  include  a  rate  effect.  In  a  dipolar  mater¬ 
ial  it  is  assumed  that,  in  addition  to  the  usual  traction  forces  that  do 
work  against  velocities,  there  are  hypertractions,  denoted  a  and  having  the 
dimensions  of  stress  times  length,  that  do  work  against  velocity  gradients. 
Decomposition  of  strain  gradients  into  elastic  and  plastic  parts  introduces 
a  new  internal  variable,  denoted  dp  in  this  paper,  which  requires  its  own 

evolutionary  constitutive  equation.  As  shown  in  Reference  9,  dipolar  plas¬ 
ticity  introduces  length  scales  that  are  characteristic  of  the  material 
itself.  In  this  paper,  it  is  supposed  that  there  is  only  one  intrinsic 
length  scale  for  elastic,  plastic,  or  viscoplastic  deformation,  which  is 
denoted  1.  In  addition  to  the  nondimensional  quantities  introduced  previ¬ 
ously,  we  require 

dp  =  hdp,  a  =  o/Ikq,  i  =  i/h  . 

The  full  set  of  equations  with  a  Von  Mises  type  flow  condition  now 
takes  the  form 

Momentum:  (s  -  £o,y),y  =  pv 
Energy:  fl  =  ke,yy  +  Srp  +£ddp 


Constitutive  Equations: 

S  =  u(v,y  -  yp)  ,  Yp  =  AS 


5  =  -  V  •  “p  ■  I  0 


Reference  Response  and  Work  Hardening: 


<=(l  +  ^-)n  ,  W  =  =  S  y  +  2nd 


Yield  condition: 

(s2  +  o2)1/2  =  (1  -  ae)(l  +  bA ( s2  +  o2)1/2)m  *  . 

In  these  equations  A  plays  the  role  of  a  plastic  multiplier  and  is  to  be 
determined  from  the  yield  condition  when  the  loading  point  (s,  o,  0)  lie3 
outside  the  static  yield  surface,  as  determined  by  the  last  equation  with 
A  set  equal  to  zero.  Extra  boundary  conditions  are  also  required.  These 
are  chosen  to  be  n(*1,t)  =  0,  so  that  the  dipolar  effect  will  effectively 
vanish  outside  the  shear  band. 


SR 
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Homogeneous  solutions  to  these  equations  are  exactly  the  same  as 
before,  but  now  the  response  to  perturbations,  computed  for  case  B  with  the 
smaller  temperature  bump,  appears  to  be  quite  different,  although  long  com¬ 
putational  run  times  have  not  yet  been  achieved.  Figure  8  shows  the  evolu¬ 
tion  of  the  central  plastic  strain  rate  for  a  nondimensional  intrinsic 

_p 

length  Z  2  10  .  The  rate  decays  rapidly  back  toward  the  average  applied 

rate,  at  least  for  short  times  after  introduction  of  the  perturbation,  so 
that  the  dipolar  effect  appears  to  stabilize  the  deformation. 


IV.  DISCUSSION  AND  CONCLUSIONS 

The  response  to  the  small  temperature  perturbation  i3  qualitatively 
very  similar  for  both  oases  A  and  B.  Figure  3  shows  that,  although  the 
stress  levels  are  different  due  to  the  different  viscous  response,  the 
shapes  of  the  homogeneous  response  curves  are  similar,  and  so  are  the 
shapes  of  the  response  to  perturbation.  For  something  on  the  order  of  5  or 
6t  average  strain  past  the  peak  of  the  homogeneous  curve  the  stress  path 
deviates  only  slightly  until  stress  collapse  sets  in  at  the  end  of  the  cal¬ 
culation.  Similar  behavior  was  reported  by  Merzer6  who  used  a  completely 
different  flow  law  for  viscoplasticity.  Thus,  it  appears  that  stress  col¬ 
lapse  is  inherent  in  the  process  of  band  formation  with  the  details  depend¬ 
ing  on  the  specific  material  model  used. 

The  calculation  with  the  larger  perturbation  in  temperature  for  case  A 
shows  that  the  details  of  stress  collapse  also  depend  strongly  on  the  mag¬ 
nitude  of  the  perturbation.  Although  case  B  has  not  been  run  with  the 
larger  temperature  disturbance,  there  is  no  reason  to  expect  that  the 
results  would  be  qualitatively  different  from  case  A. 

The  most  striking  aspect  of  the  response  patterns  in  Figure  3  is  their 
similarity  to  those  obtained  in  analyzing  imperfection  sensitivity  for  a 
bifurcation  problem,  e.g.,  see  Reference  10.  The  only  thing  missing  is 
the  bifurcation  branch  itself.  Figures  and  5  for  temperature  and  plastic 
strain  rate  reinforce  this  idea. 

Although  calculations  for  the  dipolar  case  have  not  yet  progressed  far 
enough  to  be  definitive,  it  appears  that  the  effect,  even  for  a  very  small 
inherent  material  length  scale,  will  be  highly  stabilizing.  This  would  be 
in  accordance  with  the  experience  to  date  for  elastic  systems  where,  as  in 
Reference  11,  higher  order  effects  allow  smooth  transitions  instead  of  dis¬ 
continuities. 

Unfortunately,  because  of  computational  instabilities  that  set  in,  it 
has  not  yet  been  possible  to  continue  calculations  to  the  point  where  a 
final  configuration  for  the  shear  band  appears.  The  calculations  reported 
in  Reference  2  indicated  that  a  peak  value  for  plastic  strain  rate  may 
occur,  but  reexamination  has  indicated  numerical  instability.  The  diffi¬ 
culty  seems  to  be  that  eventually  the  extreme  localization  of  the  band 
defeats  the  ability  of  the  present  numerical  method  to  resolve  the  details. 
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